Electronic structure of solid coronene: differences and commonalities to picene 
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We have obtained the first-principles electronic structure of solid coronene, which has been recently 
discovered to exhibit superconductivity with potassium doping. Since coronene, along with picene, 
the first aromatic superconductor, now provide a class of superconductors as solids of aromatic 
compounds, here we compare the two cases in examining the electronic structures. In the undoped 
coronene crystal, where the molecules are arranged in a herringbone structure with two molecules 
in a unit cell, the conduction band above an insulating gap is found to comprise four bands, which 
basically originate from the lowest two unoccupied molecular orbitals (doubly-degenerate, reflecting 
the high symmetry of the molecular shape) in an isolated molecule but the bands are entangled as 
in solid picene. The Fermi surface for a candidate of the structure of K x coronene with x = 3, for 
which superconductivity is found, comprises multiple sheets, as in doped picene but exhibiting a 
larger anisotropy with different topology. 

PACS numbers: 74.20.Pq, 74.70.Kn, 74.70.Wz 



Introduction — Discovery of superconductivity in solid 
picene doped with potassium atoms [l| is seminal and 
came as a surprise in that the picene, on top of the high- 
est T c among organic superconductors, belongs to aro- 
matic compounds, the most typical, textbook class of 
organic materials. Then a natural question to ask is: can 
other aromatic compounds become superconducting as 
well? Before examining this question, let us first briefly 
summarize the known superconductors in carbon-based 
materials. First discovery goes back to 1965, when a 
graphite-intercalation compound (GIC), KCg, was found 
to be a superconductor at 0.1 K and the highest T c 
among GIC's to date is 11.6 K in CaCeQ. Fullerene is 
another class, where potassium-doped fullerene, K3C60, 
has T c = 20 K [4], followed by Cs 2 RbC 60 with T c = 33 
K I, or 40 K in Cs 3 C 60 under 15 kbar 0. In the 
last decade, the boron-doped diamond 0] joined carbon- 
based superconductors. The first aromatic superconduc- 
tivity in the potassium-doped solid picene discovered by 
Kubozono's group[l|, with T c = 7 — 18 K in K^picene 
for x ~ 3, was unsuspected, since the (undoped) solid of 
picene (a hydrocarbon compound with five benzene rings 
connected in a zigzag) has been known to be a good in- 
sulator, as naturally expected for an aromatic molecule. 

The present authors have reported for the first time 
the first-principles electronic structure of both undoped 
and doped solid picene in the framework of the density 
functional theory (DFT) within the local density approx- 
imation (LDA) Q. We have revealed, first for the un- 
doped solid picene, where picene molecules take a her- 
ringbone structure, that the conduction band consists of 
four bands, which originate basically from the lowest two 
unoccupied molecular orbitals (LUMO and LUMO+1) 
of an isolated molecule, but the bands are entangled ( 
i.e., crossing with each other. When doped with potas- 
sium atoms, the herringbone structure is deformed, and 



the electronic wave function significantly spills from the 
organic molecules' LUMO's into potassium sites |8j. 

Now, if we come back to the question of whether and 
which other aromatic compounds can become supercon- 
ducting, recently Kubozono's group reported supercon- 
ductivity in doped coronene [9j]. Coronene, C24H42, is 
another of typical aromatic molecules, with seven ben- 
zene rings assembled in a concentric disk. Supercon- 
ductivity in K x coronene is reported to appear around 
x ~ 3 with T c up to 15 K.[Io[ Motivated by this, here we 
have obtained the electronic structure of solid coronene, 
both undoped and doped. An obvious interest is the 
differences and commonalities between the crystal and 
electronic structures of solid coronene as compared with 
those of solid picene. This is precisely the purpose of 
the present work. The obtained results are analyzed 
with maximally- localized Wannier functions (WF's) [ll|, 
in terms of which we downfold the system into a tight- 
binding model and compare with those for picene. 

We shall show that the conduction band of the un- 
doped solid coronene comprises four bands, which basi- 
cally originate from the two LUMO's (doubly-degenerate, 
reflecting the symmetry of the molecule higher than that 
of picene) in an isolated molecule, but the bands are en- 
tangled as in solid picene. The Fermi surface for a candi- 
date of the structure of K^coronene with x = 3, for which 
superconductivity is found, comprises multiple sheets, as 
in doped picene but exhibiting a more one-dimensional 
character with different anisotropy nad topology. 

The present electronic structure will serve as a ba- 
sis for discussing mechanisms of the superconductivity. 
Since some classes of organic superconductors are con- 
sidered to have an electronic mechanism, [12[ picene and 
coronene superconductors may possibly belong to them. 
There have been some theoretical studies pjj Q that 
suggest that solid picene is a strongly correlated elec- 
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tron system. As for the electron-phonon coupling, on the 
other hand, Kato et al. estimated the phonon frequen- 
cies and the electron-phonon coupling for various hydro- 
carbon molecules [15]. For picene they found that the 
electron-phonon coupling can be as large as 0.2 eV. While 
it is generally difficult to quantitatively estimate T c for a 
phonon mechanism due to the ambiguity in the Coulomb 
pseudopotential (/i*), they concluded that T c ~ 10 K 
may be expected. They also estimated the coupling for 
coronene to be about 0.1 eV. Subedi et a/. [16] recently 
calculated the phonon spectrum and electron-phonon in- 
teraction for solid picene, and they found that the calcu- 
lated electron-phonon coupling 0.1-0.2 eV is sufficiently 
strong to reproduce the experimental T c of 18 K within 
the Migdal-Eliashberg theory. 

First-principles bands - - The calculation is based 
on DFT, where LDA in the Perdew-Zunger for- 
mula is adopted for the exchange-correlation energy 
functional. [13 We use the projector- augmented wave 
(PAW) method [HI, implemented to the Quantum MAte- 
rials Simulator (QMAS) package. [Hj The pseudo Bloch 
wave functions are expanded by plane waves up to an 
energy cutoff of 40 Ry with 4x6x4 /c-points. 

As for the crystal structure, here we adopt the experi- 
mental lattice parameters of natural coronene (karpatite) 
reported by Echigo et al. [20| Natural coronene has a 
monoclinic (space group: P2\/a) structure with a = 
16.094,6 = 4.690 and c = 10.049 A, and f3 = 110.79°. 
We have also performed calculations adopting the lat- 
tice parameters of synthetic coronene [2l| , and have con- 
firmed that the results are essentially unchanged from 
those for the natural coronene. The molecular solid has 
a herringbone arrangement of molecules as depicted in 
Fig. [1] with a unit cell containing two molecules (cen- 
tered respectively at (0, 0, 0) and (1/2, 1/2, 0) as dictated 
by the symmetry). The lattice parameters are fixed at 
the experimental values, and the internal atomic posi- 
tions are optimized. The angle between the planes of the 
inequivalent molecules in the optimized geometry is 95°, 
which agrees with the measured value. It was found that 
the point-group symmetry D§h for an isolated coronene 
molecule is lowered to D2h in the crystal. 

Figure 2 displays the electronic band structure of 
coronene. We have an insulator with a band gap of 2.41 
eV. The gap, which is naturally smaller than the HOMO- 
LUMO gap, calculated to be 2.90 eV, of an isolated 
coronene molecule, is indirect, with the valence band top 
located at Y in the Brillouin zone, while the conduction 
band bottom at V. Since the crystal structure is layered, 
where the herringbone arrangement, on the a — b plane, 
of the molecules are stacked along c axis, the electronic 
structure is more dispersive along a* — 6* axis. More 
precisely, the dispersion along 6* is larger than that for 
a*, which directly reflects the distance between neigh- 
boring molecules being shortest along b. The conduction 
band, with a width of 0.40 eV, consists of four bands de- 
rived from the doubly-degenerate, e± g LUMO's. For solid 
picene by comparison, the conduction band, with a width 



0.39 eV, consists of four bands with a band gap of 2.36 
eV, while the LUMO-HOMO gap of a picene molecule 
is 2.96 eV. These are similar to solid coronene, where a 
difference is the conduction band of the solid picene orig- 
inates from LUMO and LUMO+1 of an isolated picene 
molecule Q. Considering the general tendency of LDA 
to underestimate band gaps, we expect that the actual 
band gap of coronene may be larger than the calculated 
value. An accurate prediction of optical properties and 
band widths will also require incorporation of many-body 
effects, as is done with GW by Roth et al. [22[ for picene. 
The valence band of solid coronene has a width of 0.45 
eV, and consists of four bands, which is shown to be de- 
rived from the (again doubly-degenerate) C2 U HOMO's of 
an isolated coronene molecule. 
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FIG. 1: (Color online) Crystal structure of undoped solid 
coronene viewed along (a) b axis, (b) c axis and (c) a bird's 
eye view. Solid lines delineate unit cells. 

Downfolding — We have constructed an ab initio tight- 
binding model from the maximally localized WF's PJJ of 
the four-fold conduction band. The two-fold degenerate 
e\ g LUMO's for an isolated molecule have bs g and b2 g 
symmetries in the D2h representation. Figure [3^a) dis- 
plays the two WF's localized at each of the molecules. 
While the downfolded tight-binding band (not shown) 
accurately reproduces the DFT-LDA band structure, we 
can notice that the WF is derived mainly from bs g 
orbital while w\ from b2 g . The difference in the orbital 
energy between the two WF's is 18.2 meV, which is to be 
compared with 11.3 meV for solid picene. Figure [3](b,c) 
depicts the major transfer integrals in the downfolded 
tight-binding model, which exhibit the transfers along b 
axis much stronger than the other axes, in contrast to 
those in picene (see Fig. 3(c) in Ref.Jij]), for which two 
transfers on a — b plane are close to each other. The 
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FIG. 2: (Color online) Calculated electronic band structure 
of the undoped crystalline coronene. The origin of the en- 
ergy is set to be the valence band top. The inset depicts the 
Brillouin zone with T, Z, C, Y, X, D, E and A respectively 
corresponding to (0,0,0), (0,0,1/2), (0,1/2,1/2), (0,1/2,0), 
(1/2,0,0), (1/2,0,1/2), (1/2,1/2,1/2) and (1/2,1/2,0) in 
units of (a*, b*, c*). 
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FIG. 3: (Color online) (a) displays maximally localized Wan- 
nier functions (right) wi (wh) with a lower-energy (higher- 
energy) constructed from the conduction bands of undoped 
solid coronene. Their transfer integrals t on the a — b plane 
(b) and b — c plane (c) are also shown, where the red (blue) 
arrows are for \t\ > 30 meV (10 < \t\ < 30 meV). 



anisotropy of the transfers in coronene thus comes from 
the relative distances between the WF's. 

Doped solid coronene — Let us finally discuss the 
doped coronene. If we naively adopt a rigid-band pic- 
ture for a doping level of x = 3 for which superconduc- 
tivity is observed, the Fermi surface (not shown) con- 
sists of multiple surfaces that comprise one-dimensional 
(planar) surfaces and a two-dimensional (cylindrical) 
one, along with a more three-dimensional one. How- 
ever, since we have shown for the doped picene [8j that 
the rigid-band picture is broken in this material for 
two- fold reasons (i.e., a distorted herringbone structure 
upon doping along with a spilling of the molecular wave 
function over to potassium sites), we should have a 
look at the electronic structure of the coronene actu- 
ally doped with K. Since the structure of the doped 
system has not been experimentally obtained, we have 
tried a structural optimization of Kscoronene. Since 
this in itself is an important but vast task, here we 
only show a candidate structure. We started from 
a plausible geometry with one potassium atom just 
above the molecule and two in the interstitial region. 
The initial coordinates of the six potassium atoms we 
adopted are (1/2, 0, 0), (0, 1/2, 0), (±1/3, 1/2, ±1/3) and 
(±5/6, 0, ±1/3), which preserve the monoclinic symme- 
try. Interestingly, the optimization, with the lattice pa- 
rameters fixed at those for the pristine structure, resulted 
in a significant rearrangement of the herringbone struc- 
ture with the monoclinic symmetry lowered, along with a 
strong deformation of each molecule, as seen in Figs. 2{a) 
and (b). The band structure of the Kscoronene is signif- 
icantly more dispersive than the undoped one, where the 
LUMO-derived band group is fused with the upper band 
group, resulting in a much wider band group, as shown 
in Fig. IHc). We can see that the resultant Fermi sur- 
face is an intriguing composite of one-dimensional sur- 
faces (i.e., multiple pairs of planar surfaces) with dif- 
ferent anisotropics. The doped picene also has a co- 
existence of multiple Fermi surfaces but the doped 
coronene has topology of the surface different from those 
in doped picene, which comes from different (and more 
one-dimensional) tight-binding structure (Fig. [3^b,c)). 
The anisotropic, multiple Fermi surface should give a ba- 
sis for examining superconductivity. In a phonon mecha- 
nism the problem becomes the coupling between the elec- 
trons on such a Fermi surface and the molecular phonons. 
The situation gives an interesting possibility for electron 
mechanisms as well, since the nesting between discon- 
nected Fermi surfaces can give rise to a unique opportu- 
nity for an electron mechanism, especially in multiband 
cases. [23] For an accurate description of the structure 
more elaborate and exhaustive structural optimizations 
will be needed, since the large degrees of freedom on the 
dopant positions will render the energy surface many lo- 
cal minima. Such a systematic examination is, however, 
beyond the scope of the present work, and will be re- 
ported in future. 

Finally, a few words about the "aromaticity" of 



molecules. Its simplest definition is in terms of the 
"Clar sextets" (resonating benzene rings) [24j on a given 
molecular structure. When a molecule is fully benzenoid 
(where the sextets exhaust all the double bonds) the wave 
function tends to be localized on the sextets. In this pic- 
ture picene is non-fully-benzenoid, while coronene, also 
non-fully-benzenoid, has multiple Clar structures. (2l) 
The relation of these quantum chemical properties with 
band structures is another of interesting future problems. 

In summary, we have studied the electronic structure 
of solid coronene by means of first-principles calculation. 
The conduction band is found to comprise four bands, 
which basically originate from the lowest two unoccu- 
pied molecular orbitals (doubly-degenerate, reflecting the 
molecular symmetry) in an isolated molecule, but the 
bands are entangled. The maximally localized WF's are 
used to derive a downfolded tight-binding Hamiltonian, 
where the major transfer integrals are found to signifi- 
cantly differ from those in solid picene. The Fermi sur- 
face for a candidate of the structure of K^coronene with 
x = 3, for which superconductivity is found, comprises 
multiple sheets, as in doped picene but with different 
topology of the surface. Their relevance to superconduc- 
tivity is an interesting future problem. 
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FIG. 4: (Color online) Crystal structure of K3 coronene 
viewed along (a) b or (b) c axis. Larger balls represent K 
atoms, (c) shows electronic band structures of the K3 coronene 
(circles) and undoped coronene (curves). The origin of en- 
ergy for the Kacoronene is set to its Fermi level while that for 
undoped coronene set for x = 3 in a rigid-band shift for com- 
parison. The Fermi surface for the K3 coronene is displayed 
in (d). 
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